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ABSTRACT 

The  stability  of  pseudospectral-Chebyshev  methods  is  demonstrated 
for  parabolic  and  hyperbolic  problems  with  variable  coefficients.  The 
choice  of  collocation  points  is  discussed.  Numerical  examples  are  given 
for  the  case  of  variable-coefficient  hyperbolic  equations. 
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1. 


Introduction 


The  purpose  of  this  paper  is  to  analyze  spectral  Chebyshev  collocation 
(otherwise  known  as  pseudospectral)  methods  for  hyperbolic  and  parabolic 
problems.  We  shall  show  that  these  methods  converge  at  a  rate  that  is 
faster  than  that  of  finite  differences.  The  analysis  is  based  upon  results 
presented  in  [l].  This  reference  outlines  the  general  theory  of  convergence 
of  spectral  methods  and  proves  that  if  a  spectral  method  is  algebraically 
stable  in  some  norm,  then  the  method  is  strongly  stable  in  an  algebraically 
equivalent  new  norm.  If  in  addition  the  method  is  consistent  by  virtue  of 
its  truncation  error  tending  to  zero  in  this  new  norm,  then  convergence  is  implied. 

The  application  of  this  theory  to  hyperbolic  and  parabolic  problems 
had  been  discussed  in  [l]  mainly  for  constant  coefficient  hyperbolic  and 
parabolic  problems  and,  in  the  case  of  Chebyshev  methods,  mainly  for  the 
Galerkin  and  Tau  methods.  In  this  paper  we  discuss  the  collocation 
methods  and  prove  stability  for  the  variable  coefficient  case.  The  new 
idea  that  enables  us  to  establish  stability  for  collocation  methods  is  a 
new  quadrature  formula  for  Gauss-type  integration.  We  use  the  positive  weiehts 
given  by  this  formula  as  the  new  norm  and  prove  energy  conservation  in  this 
norm.  Using  the  same  technique  a  new  proof  is  presented  for  variable 
coefficient  hyperbolic  and  parabolic  problems  when  solved  by  spectral- 
Chebyshev  methods  using  Tau  methods.  These  proofs  are  more  general  than 
those  in  [l]  in  the  sense  that  they  include  the  variable  coefficient  case. 


Section  1 

A  numerical  solution  of  the  problem 

ut  ■  Lu 


(1.1) 


where  u  £  H,  H  is  an  Hilbert  space  and  L  is  an  infinite  dimensional  operator » 
consists  of  two  steps.  The  first  is  to  choose  a  finite  dimensional  sub¬ 
space  of  H,  say  3^,  and  the  second  is  to  choose  a  projection  operation 
Pjj  :  H  Bn«  The  approximation  to  (1.1)  becomes 


p  L  P 
N  N 


N 


u  €  B 
N  N 
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which  may  be  solved  on  a  computer.  Spectral  Chebyshev  methods  are 
defined  by  choosing  B^  as  the  N-dimenslonal  space  spanned  by  polynomials 

of  degree  N+k-1  that  satisfy  boundary  conditions. 

There  are  three  ways  which  have  been  used  to  choose  the  operator  , 

namely  Galerkin,  Tau  and  collocation. 

In  the  Galerkin  method  for  homogeneous  boundary  conditions  we  choose 

4>  a  -  1, ...  ,N  as  the  basis  of  B„  and  solve 
n  N 


(  \ 

\  L  V  ♦n  /  *  °  n  "  1 . N 

(1.3) 

N 

“N  "  E  an  ’ 


For  the  Tau  method  we  choose  Co  be  a  set  of  orthogonal  functions 

such  that  (0  ,0  )  ■  5 _  and  expand 

n  m  tun 

N+k 

"s  "  L  a„ 

n«l 


where  k  is  the  number  of  boundary  conditions.  Then  set 
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(£-‘v  J 


n  *  1, . . .  ,N 


Cl. 4) 


The  condition  6  provides  the  other  k  equations. 

In  the  collocation  method  we  set 


N 

I 

n=*l 


a„<|> 
n  n 


and  require 

3u 

^  =  0  for  xj  j  =  1 . N  •  (i- S) 

It  had  been  observed  by  Orszag  [  1  ]  and  Kreiss  and  Oliger  [  2  ]  that  the 
collocation  method  can  be  carried  out  efficiently  in  the  physical  space  in 
contrast  to  the  Galerkin  and  Tau  methods  which  must  be  solved  in  the 
transform  space.  This  fact  enables  one  to  use  the  collocation  method 
efficiently  for  nonlinear  equations.  We  refer  the  reader  to  [l]  for  further 
discussion  of  this  fact. 

In  the  next  sections  we  will  illustrate  the  above  procedure  applied  to 
parabolic  and  hyperbolic  equations. 


Section  2 

Consider  the  equation 

-1  <  x  <  1 
0  <  6  <  S(x) 


S(x)  u 


XX 


u(±l,t) 


(2.1) 
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In  chc  Galerkin-Chebyshev  method  we  choose 


■ 

T 

-  T. 

n  even 

n 

n 

0 

* 

T 

-  Ti 

n  odd 

n 

T  (x)  »  -os(n  cos  ”^x) . 

n 

1 

\'o  oxnand 


“s  ■  Z  a»  +n(x) 
a»0 


so  that  u^(±I,t) 


0  and  set 


(2.2) 


It  is  readily  seen  that  for  nonconstant  S(x) ,  it  is  difficult  to  solve  the 

equations  for  (2.3).  Orszag  has  found  some  efficient  transform  methods  to 
evaluate 


In  general,  however,  solving  (2.3)  for  the  coefficients  (an)  1*  Mum  consuming 
In  the  Tau  method  we  sec 


N+2 

s-  5. 


n-0 


T  (x) 
n 


and  require 


S(x)  — y-\  T  dx 
3x2  /  n 


0 


n  ■  0, 


-1 


We  face  the  same  complications  for  getting  the  coefficients  as  we  had  for  the 
Galerkin  method. 

In  the  collocation  method  we  set 


N 


•  ■  UH  ‘  £.  a»  *n(x) 
n=0 


where  the  4,n's  are  defined  in  (2.2).  Then  we  demand 


3u 


N  ,  N  '  “n 

—  -  S(x)  71"  =  °  3t  X  =  Xj 


9x 


j  =  0, . . .N 


(2.5) 


for  some  points  x ^ .  If  the  x^  are  chosen  to  be  cos  so  that  the 


boundary  values  are  included,  there  is  an  efficient  way  to  solve  (2.5),  by 
taking  advantage  of  the  orthogonality  of  the  trigonometric  functions.  Set 


Vx)  *  £  Vn(x)  *  W  =  S  anTn(V  ’  °-j-N 

n=0  n=U 
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Then 


N 

2  V'  1  ,  .  Trnk 

a  *  —  2-i  —  u.t(x.  )  cos  — — 

a  Me  ,n  c,  N  k  N 

n  k»0  k 


=0  ■  CN  '  2 

ck  -  1  l<k<N-l 


(2. 


Using  Che  properties  of  Chebyshev  polynomials  we  may  set 


3  u  (x. )  N 
'  /  3  -  E  b„  T  or.) 
3x2  nto  n  n  J 


where  the  coefficients  may  be  found  from 


N 

b  «  T 

in 


2  2 

p(p -n)a 


p  *  n+2 
p+n  even 


Then  we  go  back  to  the  physical  space  and  solve 
3u)i  3%. 

T^kx,)  -  S(x.)  —A  (*,)  j  ■  1 . N-l 

ac  j  J  3xz  3 

W  “  W  *  0  • 


This  procedure  is  very  efficient  and  may  be  generalized  without  any 
problem  to  nonlinear  equations.  In  practice  we  would  use  the  Chebyshev 
polynomials  to  interpolate  u  spatially  and  then  to  evaluate  the  spatial 
derivative  at  the  desired  points  x^ .  Finally  the  solution  would  be 
advanced  in  time  using  the  original  nonlinear  equation  to  find  the  time 
derivative  at  the  points  Xj  in  the  physical  space. 

In  order  to  prove  convergence  we  need  the  following  two  results. 


5a) 


-6- 


Lemma  1 . 


Let  u  satisfy  u(±l)  -  0  and  have  a  continuous  first  derivative,  then 


dx  <  0 


(2.6) 


For  the  proof  we  refer  the  reader  to  [1,  p.  82], 


Lemma  2. 

Let  x.  =  cos  „  j  =  0, 
J  N 


,N.  Then  there  exist 


w^.  >  0,  j=0, . 


.N 


such  that 


N 


I 

j=0 


f(x  ) 


w . 
J 


(2.7) 


w . 


j=0 


7T 


for  any  f(x)  which  is  a  polynomial  of  degree  at  most  2N-1. 

Remark:  The  formula  (2.7)  is  a  generalization  of  the  well  known  integra¬ 
tion  rule  of  Gauss  type.  Note  that  depends  on  N.  There  are  two  major 

differences  between  (2.7)  and  the  usual  Gauss  integration  formula  for  the 

2  -x 

weight  function  w(x)  =  (1-x  )  2.  The  first  difference  is  that  the  Gauss 

formulas  are  of  open  type,  that  is  the  boundary  points  are  not  included, 

whereas  in  (2.7)  they  are  included.  The  second  is  that  the  interior  points 

are  not  the  zeroes  of  the  orthogonal  polynomials  with  respect  to 
2  _i 

w(x)  =  (1-x  )  namely  the  Chebyshev  polynomial  of  the  first  kind,  but 
rather  that  they  are  the  zeroes  of  the  Chebyshev  polynomial  of  the  second 
kind  (which  are  orthogonal  with  respect  to  the  weight  function 

2  h 

w(x)  =  (1-x  )  ).  It  is  interesting  to  note  that  Lemma  2  implies  the  Gauss 

2  4 

integration  formula  for  the  weight  function  w(x)  =  (1-x  )  ,  i.e.,  that 


the  formula 


N-l 

j-1  J 


(2.8) 


is  correct  for  any  2N-3  degree  polynomial.  In  fact  if  (2.7)  is  correct 
2 

set  f(x)  -  (1-x  )g(x)  in  (2.7)  to  get 


since  Xq  -  1,  x^  *  -1,  and  the  fact  that  if 

polynomial  then  the  degree  of  f  is  2N-1. 

2 

established  with  w^  -  (1-x^)  w^  . 


N-l 

-  2Z  (l-x.)g(x.)  w 
j-1  J  J  J 

g(x)  is  a  2N-3r<*  degree 
Equation  (2.8)  is  now 


Proof: 

Wo  first  note  that  (2.7)  can  be  made  exact  for  any  NCtl  degree  polynomial 
by  putting  f(x)  »  xQ,  n‘--0,...,N  and  solv-1  -g  for  the  N+l  unknowns, 
w_.  since  the  Vandermonde  matrix  is  aonsingular. 

Let  f(x)  be  now  a  polynomial  of  degree  2N-1.  Then  there  are  g(x) 
of  degree  N-2  and  v(x)  of  degree  N  such  that 

f(x)  -  (l-x2)YJ{_1(x)g(x)  +  v(x)  (2.9) 


where  Y^_^  is  the  Chebyshev  polynomial  of  the  second  kind,  i.e., 


Since 


sinCNcqs~Vx) 


N-l 


a-x 


(2.10) 


V 

kN-l 


(x.) 


mi  = 

zi 


j-1, . . .N-l 
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we  conclude  that 


f(x.)  =  v(x.) 
J  J 


j=0, . . . ,N 


(2.11) 


1  1  _  1 

f  -  -X—  dx  =  f  /l-x2  YN_(x)g(x)dx  +  f  AAL.  dx  #  (2  12) 

J  J  J.  .X  _„2 


-1  ^1- 


The  first  term  in  the  right  hand  side  of  (2.12)  vanishes  since  is 

orthogonal  to  any  polynomial  of  degree  less  than  N-l  and  the  degree  of 
g(x)  is  N-2.  Since  v(x)  is  a  polynomial  of  degree  N  (2.7)  is  exact 
for  it  and  therefore 


_1  /l-x 

and  by  (2.11)  we  conclude  that 


=  L  v(xi>  wi 
j-0  J  3 


f(X)  dx  = 

Z*1  v(x) 

A-x2 

l  £? 

This  proves  (2.7)  except  for  the  fact  that  w  >  0.  To  prove  that  we  define 


fz(x)  =  ( l-x2)  Y2_1  (x)  /  [  Y^j_1(x^)(x-x^)  ]2  l  =  1,...,N-1  . 


The  degree  of  f^(x)  is  2N-2;  moreover,  =  <5£j^1_X  ^  and 

f.(:<)  >0,  -1  <  x<  1.  Therefore  (2.7)  is  exact  and  yields 


-9- 


r  Vx)  iL  2 

/  -=.  -  Z  W  »j  -  <1-* )  «t 

'  L  /l-x“  j*0 


(2.13) 


Equation  (2.13)  shows  that  w.  >  0,  l  »  1,...,N-1.  Define 


now 


and 


fg(x)  »  (1+X) 

N-l 

(a) 

2  -  0 

fN(x)  »  (1-x) 

(b) 

(2.14) 


W  ■  % 

and  therefore 


W  '  5Sj 


w  *  /  ■  1  ■  dx  >  0 

4  /£? 

wv  *  /  dx  >  0 

4  £7 


and  this  concludes  the  proof. 


(a) 


(b) 


(2.15) 


We  are  now  ready  to  prove  the  stability  of  the  Chebyshev  collocation  method 
for  Che  heat  equation. 


Theorem:  (Stability) 

Let  u^  be  Che  Chebyshev  collocation  a, proximation  (2.5)  to  the  heat 
equation  (2.1).  Then 


N  u^(x.t)  N  ,  w. 

Elfrh-i  E  ^(xr0)l( 

J-0  ^  V  J  j-0  J 


x.) 


(2.16) 
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Proof : 


Since 


3u 

_ H 


.2 

9  UN 

(x  ,t)  =  S(x. )  — r—  (x.)  j=l,...,N-l  and 


and 


We  get 


3t  '  j 

W  =  W  "  0 


j  ax2  j 


N  3u  w .  N 

2  W  sr  ‘"j1  lory 


a2 

2  W  ~T  (xj}  wj 

j=0  J  3x  J  j 


(2.17) 


-2 

9  ^ 

By  Lemma  2.2  and  the  fact  that  the  degree  of  — =—  is  2N-2  we  get 

3x 


..  22  1  -2 
N  -  3  uv  r  ^  3  Sj 

j?o  VV^r(V“j-/ 

J  _]_  /1-x 


dx  <  0 


The  last  inequality  follows  from  Lemma  2,  since  u^(±l,t)  =  0.  Therefore, 


“n'V 


W  . 

J- 


s(xj) 


<  0 


and  (2.16)  follows. 

The  next  step  for  showing  convergence  is  to  show  that  the  truncation 
error  tends  to  zero  as  N  ^  tends  to  infinity.  In  view  of  the  discussion 
in  [1,  p.  48]  the  truncation  error  is  given  by 


PH  SW  7-2  W00  7-2 

3x  3x 


(2.18) 
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where  u  £  C*  is  the  solution  to  (2.1),  Pjjf(x)  is  a  polynomial  of  degree 
N  that  interpolates  the  function  f(x)  at  the  points  x^  and 


N  g  (x.) 


z 

U-o 


1 


s(x4)  Wj 


cos 


j-0, . . . ,N 


Theorem:  (Consistency) 

Let  u,  and  j| 


be  defined  as  above  then 


2  2 
PN  S(x)  77  PN~  S(x)  77  u 

9x  3x 


(2.19) 


for  any  positive  r. 


Proof: 

From  (2.5a)  we  can  express  PNu  by 

N  a 

V  ■  £  r  V*> 

n-0  n 


where 


,  N  u(x  ) 

-  y  — i- 

j-0  cj 


Vxj> 


On  the  ocher  hand 


where 


.  r  u(x)T  (x) 

,  -  ^  / - 2 d* 

"  \J  /T? 


It  is  well  known  that  a  ■  0 

n 


for  any  p.  Moreoever  a  can  be 

n 


expressed  in  terms  of  the  an's  by  the  formula 


.  N  T,  (x.)T  (x.)  ,  N  NT,  (x.)T  (x.) 

„  aIa  k  J  n  i'  ,  2  y*  a  V"  k  J  n  J 

-  N  a0  i.  ”  ^ 

J=0  J 


EakE 

k=l  j=0 


=  an  +  a2N-n  +  a4N-n  + 


and  therefore 


S(x)  2  PNU 

dX 


S(x)  fl  u  =  S(x)  £  (an  ~  an)Tn 
3x  n=0 


-  s<*>  £  v; 

n=N 


define 


then 


g(x)  =  S(x)  £  a  T" 
n=N 


j=0 


00 

,2 

r «, 

E 

n=N 

a  T"  (x.) 
n  n  j 

J 

J* 

H 

O 

1 

Nr 

Moreover  since  a  -a 
n  n 


by  (2.20)  we  get 


(2.20) 


(2.21) 


and  since 


II  S(x)  2—  PNu  -  S(x)  1—  u 
3x  3x 

*  1  (2.19)  is  proven. 


Section  3 

In  this  section  we  would  like  to  treat  the  hyperbolic  equation 


u£  -  S(x)ux  S(x)  >0  | x |  £  1 

u(l,t)  -  0  . 


(3.1) 


We  concentrate  upon  the  collocation  method.  There  are  currently  two 

ways  of  performing  the  collocation  method.  The  first  one  is 

irk 

to  collocate  at  the  point  x^^  -  cos  ~  ,  k-1 . N  and  to  use  the  boundary 

condition  for  xQ  -  1.  This  means  that  we  collocate  at  N-l  points  in  the 

interior  of  the  domain  and  also  at  the  outflow  boundary;  we  do  not  collocate 

at  x  »  1  since  a  boundary  condition  is  imposed  at  this  point.  The  other  way 

TTk 

is  to  collocate  at  the  points  x^  ■  cos  —  ,  k-1,..., N-l  and  to  use  the 

boundary  condition  at  x^  «  1.  This  amounts  to  using  N-l  interior  points 
for  collocation  and  to  impose  a  boundary  condition  at  the  inflow.  The  out¬ 
flow  boundary  is  not  treated  at  all.  We  would  now  like  to  show  how  to  carry 
out  these  two  methods  effectively. 

In  order  to  carry  out  the  first  one  we  expand 


We) 


N 

Z 

k-Q 


•k  w 


n  *  0, . . .  ,N 
im 

x  -  cos  — 
n  N 


and  solve  for  a^ 


(3.2) 
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(3.3) 


1^1 

\  m  —  E  —  ^(x  t> 

ck  j-0  cj  N  J 


cos 


TTjk 


cn"  2 


0^8.#K 


Equation  (3.3)  is  evaluated  by  using  the  Fast  Fourier  Transform  (FFT)  method. 
Now 

3un  N 

ar  (VC)  *  E  bk  Tk(xn)  .  (3.4) 

k»0 

where 

1  N 

blc  ■  ~  E  2Pap  •  (3.5) 

Ck  k-p+1  p 

k+p  odd 

The  evaluation  of  the  right  hand  side  of  equation  (3.4)  is  carried  out 
using  FFT.  Then  equation  (3.5)  is  solved  for  the  b^/s  with  0(N) 
operations, that  is  a  simple  recursive  formula  is  used 

'  bN  ‘  0  Vl  *  28an 


and 


bk+2  "  bk 


2  (k+l)a, 


'k+2 


k+1 


Then  we  solve  in  the  physical  space. 


3uN  3uN 

(Xj , t)  -  S(x^)  (Xj.t)  j  -  1.....N 


uN(l,t)  -  0  . 


(3.6) 


A  very  efficient  time  marching  technique  which  is  explicit  and  uncondi¬ 
tionally  stable  had  been  developed  in  [3]  and  can  be  used  for  the  solution 
of  (3.6). 
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The  second  way  of  collocacion  is  carried  out  as  follows.  Set 


N-l 


v  (»n,t)  -  £  d  T  (x„) 
k-0 


0 . N-l 

TO 


(3.7) 


xn  "  C0S  N 


It  can  be  shown  that  d^  can  be  expressed  In  terms  of  a^  derived 
in  (3.3).  In  fact 


\  ■  \  *  ('l)N‘l2aN  ^ 


(3.8) 


Equation  (3.8)  is  derived  as  follows 


W0 


N-l 

E  %  W  *  V»(,n) 

k-0  ' 


N-l 


E  \  VV  *  <-uS'l2aN 

k-0 


’N  ,  .,k  ,  ..N 

£  T.  (X  )  T  (x  ) 

in  c.  *  n  cxi  N  n 

k-0  k  N 


N-L  „  .  N-l 

£  a  T  (x  )  +  (-l)N_i2a  £ 

k-0  r 


,  nk  N-l 

T.  (x  )  «  £  «v  Vx„> 


,  _  c,  k  n 
k-0  k 


k-0 


for 

Now 


a  -  0 , . . . ,  N- 1  . 


3vN 

T^x  ,t) 
3x  n 


N-l 


k-0 


\  W 


where 


f  E  2”% 

ck  k-p+1  ? 

k+p  odd 


(3.9) 


and  we  solve 


3vN(x.,t) 

3t 


3vN(x  ,t) 

S(x  )  -r - — ^ - 

n  3x 


,N-1 


(3.10) 


UN(x0’t)  *  °  * 


Observe  that  In  the  second  way  of  collocation  (3.7)  -  (3.10)  is  a 

polynomial  of  degree  N-l,  whereas  in  (3.2)  -  (3.6)  it  is  a  polynomial  of 


degree  N.  The  similarity  between  these  two  different  methods  can  be  seen 

3un  3ujj 

in  the  case  where  S(x)  *  1.  Since  -j—  is  a  polynomial  of  degree  N 


that  vanishes  at  x 


-1  and  at  the  series  of  T^( x )  we  get 


3u 

Jt 

and  by  the  same  argument 

3v 

3t 


N 


N 


3um 

it  ♦ 


3x 


+  x2  T^(x) 


(a) 


(b) 


(3.11) 


It  is  interesting  to  note  that  for  the  Tau  method  one  gets  the  error  equation 


^RN 

3x~  +  T3  TN(x)’ 


(3.12) 


where  R„  is  the  Tau  approximation  to  u.  It  seems  that  the  Tau  method 
N 

can  be  viewed,  in  the  case  of  the  constant  coefficient  problem  (3.12)  as 
a  collocation  method  based  on  the  collocation  points 


cos 


1  2k-l 

2  N 


k«l,...,N 


(3.12a) 


This  observation  suggests  a  convenient  way  of  using  the  Tau  method  for  the 
variable  coefficient  case  as  well,  namely  set 


3Rm  *K 

JT  "  S(x)  3x  “  °  x  "  zk  k  ■  1, . .  •  ,N  .  (3.13) 

This  method  reduces  in  Che  constant  coefficient  case  Co  Che  Tau  method. 

In  order  Co  escablish  stability  for  che  collocation  method  described  in 
(3.7)  -  (3.10)  we  need  the  following  Lemma: 


Leiwaa 


Let  x^  -  cos  k  -  0, . . .  ,N-1  then  the  quadrature  formula 


//iS  £<:0d* 

-1 


N-l 


E  £(V 

k-0 


(3.14) 


where 


_ 1 _ 

‘Vi2 


(x-x,.)2 


(1-x) 


dx  >  0 


dx  > 


0 


is  correct  for  every  polynomial  of  degree  2N-2  or  less. 


Proof. 

Let  f(x)  be  a  polynomial  of  degree  2N-2.  Set  g(x)  •  (l+x)f(x).  Since 
g(x)  is  a  polynomial  of  degree  2N-1,  formula  (2.7)  is  exact. 


I 

f  azvw  dx  - 

-1  /  1-x* 


N 


L 

j-0 


(1+Xj)w^f (Xj) 


N-L 

£  (l+x,)w  f(x). 

j-0  J  J 


(3.15) 


_i  a_ 


Equation  (3.15)  implies  (3.14)  and  can  be  derived  by  a  standard 

argument. 

Now  let  be  the  collocation  approximation  to  u  gotten  by  (3.7)  -  (3.10). 

Then 


3t 


c(x) 


3x 


X  =  X 


n  -  1, . . . ,N-1 


V  (x  )w 
N  n  n 


Multiplying  by  —  we  get  from  (3.16) 


N-l  3v  w  N-l  .  3v 

vM(x  )  (x  )  •  V v~\  =  2  vnj (x  )  r —  (x  )  w 

",  N  n  dt  n  c(x  )  „  N  n  3x  n  n 

n-u  n  n=0 


and  by  (3.14) 


1  2 


f  /T+X  8UW  A  <  1  /  1+X  v2  I  1  i  /* _ ^ 

V/Tx  VN3Tdx^  2/rxVN  '-12  7 


-1  (1 


-x)/l-> 


(3.16) 


dx  .  (3.17) 


The  boundary  term  in  the  right  hand  side  of  (3.17)  vanishes  since 

v,(l)  =  0  and  v„  is  a  polynomial  and  therefore 
N  N 


.  N-l  ,  w 

3T  s  V*„>  Too  -  0 

n=0  n 


(3.18) 


From  Che  definition  of  w^  it  follows  that  the  norm  described  by  the 
w 

weights  j  is  algebraically  equivalent  to  the  norm  in  which  we  have 

n 

consistency,  therefore  algebraic  stability  is  proved.  The  same  idea  can. 
be  utilized  in  showing  the  stability  of  the  Tau  method.  In  fact  from  (3.12) 
it  is  evident  chat 


3RN 


J?i(1+Zj)  l-*j  R^Uj)  3e 


(2j) 


(3.19) 


where  y^  are  the  weights  in  the  Gauss-Chebyshev  integration.  From  (3.17) 
it  follows  that 


which  proves  algebraic  stability. 

The  stability  of  the  collocation  method  described  by  (3.3)  -  (3.6)  follows 
immediately  from  that  described  in  (3.7)  -  (3.10).  It  can  be  seen  from 
the  relation  (3.8).  In  fact,  setting 


at 


3x 


(1+x) 


one  gets  (3.11) (a)  from  (3.11) (b).  This  completes  the  discussion  of 
collocation  method  for  scalar  equations.  We  refer  the  reader  to  [4]  in 
which  proper  ways  of  implementing  spectral  methods  for  systems  is  discussed. 

Section  4 

The  proofs  presented  in  the  last  section  were  confined  to  the  case  in 
woi.'h  c(x)  does  not  change  sign.  This  might  be  a  weakness  of  the  theory  rather 
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than  that  of  the  collocation  method.  Numerical  experiments  using  the  pseudo- 
spectral  methods  have  indicated  that  there  is.  no  instability,  that  is  they 
show  the  solution  does  not  grow  with  N  even  when  S(x)  changes  sign. 

There  might  be  problems  owing  to  growth  in  time  of  the  solution  or  to  the 
existence  of  a  stationary  characteristic  in  the  neighborhood  of  either  the 
boundary  or  some  interior  point.  But  these  problems  seem  to  occur  because 
of  lack  of  spatial  resolution  and  not  because  of  stability.  In  order  to 
illustrate  this  fact  let  us  consider  two  equations: 

Ixl  <  1 


U_  =  -XU 
t  X 


(4.1) 


u(x,0)  =  f(x) 


and 


U.  =  XU 
t  X 


u(x,0)  =  f(x) 
u(±l,t)  =  0 


X  <  1 


(4.2) 


We  attempt  to  solve  (4.1)  and  (4.2)  by  the  Chebyshev  collocation  method. 
According  to  the  popular  belief,  there  should  be  instabilities  in  the  solu¬ 
tion  since  x  changes  sign  in  the  domain.  However,  as  indicated  in  Table 
I  below  no  such  instabilities  were  found.  As  a  matter  of  fact  we  can  prove 


theorem : 

The  Chebyshev  collocation  method  for  (4.1)  is  stable. 

Proof: 

Let  u^  be  the  Chebyshev  approximation  to  u  gotten  by  the  collocation 
method.  Then 


at 


3x 


cos 


JLi 


j«0,...,N  (a)  (4.3) 
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as  in  (3.2).  Or 


IS  + 

IS 

m  n 

3t  + 

X 

3x 

m  u 

as  in  (3.7). 

Since 

“n 

is 

xch 

polynomial. 

we  get 

3t 

where  w  is 

either  . 

"b 

or 

VN‘ 

x  »  cos 


22 


j  ■  1 , .  .  •  ,  N- 1  (b)  (4.3) 


-  x 


2s 

3x 


-1  <  x  <  1 


(4.4) 


We  now  refer  the  reader  to  ([1],  p.  85-87)  for  the  proof  that  (4.4) 
implies  stability. 


Theorem : 

The  Chebyshev  collocation  approximation  for  (4.2)  is  stable. 
Proof: 

Now  a.,  satisfies 
N 


3u, 


3x 


j*l, . . . ,N-1 


(4.5) 


Ujj(il.t)  *  0  . 

Note  that  (4.5)  has  boundary  conditions  in  contrast  to  (4.4). 
Let  Wj  be  defined  in  (2.7).  From  (4.5)  it  is  clear  that 


(V  "j 


N  Xj 


ou 


s  .  .  I  3us 


1+x.  W  IT  -(  3 T  (‘U 

j»0  j  '■ 


■>  2 


N 


(1+x)* 1-x 


(1+X)i‘(l-X)' 


(l-x+x")  £  0 


and  that  proves  stability. 

It  goes  without  saying  that  the  proofs  of  the  last  two  theorems  can 

C  (  X  ) 

be  extended  to  all  the  functions  c(x)  such  that  -  is  of  constant 

x 

sign.  We  conjecture  that  it  is  true  for  any  c(x). 

In  Table  I  we  show  the  results  of  applying  the  Chebyshev-pseudospectral 
method  to  four  equations. 

The  first  equation  is 

uc  =  (l+x)ax  |x|  £  1  (a)  (4.6) 

u(x,0)  =  sin  ttx 
u(l,t)  =  sin(2eC-l)ir 

The  solution  to  this  .problem  is 

u(x,t)  =  sin  tt  [  (l+x)eC-l]  .  (b)  (4.6) 

This  problem  has  a  characteristic  boundary  at  x  =  -1,  moreover  for 
large  t  the  solution  has  a  large  variation  in  the  neighborhood  of 
x  *  -1. 

The  second  problem  is 

ut  =  (l-x)ux  (a)  (4.7) 

u(x,0)  =  sin  nx 

The  solution  is  given  by 

u(x,t)  =  sin  rr[  l-(l-x)e  C]  .  (b)  (4.7) 
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The  line  x  ■  1  is  a  characteristic  boundary  but  in  contrast  to  the  equa¬ 
tion  (4.6)  the  neighboring  characteristics  point  from  the  boundary  towards 
the  domain.  The  third  is 


“c  ■  x  (a)  (4.8) 

u(x,0)  =  sin  fix 
u(l,t)  »  sin(TTeC) 
u(-l,t)  -  -sin(7reC) 

The  solution  is 

u(x,t)  =  sin(7rxet)  (b)  (4.8) 

And  the  fourth  is 


u(x,0) 


-x  u 

x 

sin  itx 


(a)  (4,9) 


where 


u(x,t)  *  sin(TTxe  C)  (b) 

All  these  problems  were  solved  by  Chebyshev  pseudospectral  methods  with 
che  modified  Euler  time  marching  techniques.  With  the  time  step  At  *  1/Ax 
we  advance  from  the  time  0  to  the  time  t  *  2.  Note  that  since 


min 


X"  "  *Jf-l 


N 


0 

iV 

then  Ax  .  ^  0 

'  1  ' 

l}*2; 

mm 

0 

In  Table  I  we  show  the  Chebyshev  errors  of  che  solution  of  the 


problems  (4.6)  -  (4.9).  It  is  clear  that  the  Chebyshev  collocation  method 
was  stable  for  all  these  problems  and  has  the  same  rate  of  convergence. 


However,  the  errors  for  problems  (4.6)  and  (4.8)  were  much  larger  than 
those  of  (4.7)  and  (4.9).  In  fact  taking  64  modes  in  the  solution  of  (4.6) 
and  (4.8)  produce  the  same  error  that  17  modes  produce  for  (4.6)  and  (4.8). 
This  is  a  problem  of  accuracy  and  not  of  stability.  The  question  now  is  do 
we  retain  spectral  accuracy?  To  answer  this  question  we  ran  the  problem 
(4.8)  with  smaller  and  smaller  time  steps  until  the  results  were  not 
changed  which  means  that  we  get  the  space  accuracy.  For  17  modes  we  got 
an  L2  error  of  1.16.10  \  whereas  for  33  modes  we  got  an  error  of  6.10 
This  indicates  the  fact  that  the  order  of  accuracy  in  space  is  indeed  better 
than  any  algebraic  order. 


Conclusion 

It  has  been  shown  in  this  paper  that  the  pseudospectral-Chebyshev 
methods  are  convergent  in  variable  coefficient  parabolic  problems  and  in 
some  cases  to  hyperbolic  problems.  The  analysis  shows  that  the  rate  of 
convergence  is  greater  for  finite  difference  methods  or  the  finite  element 
method.  It  seems  that  for  a  single  first  order  hyperbolic  equation  the 
method  remains  stable  even  when  the  coefficient  changes  sign,  though  in 
this  case  care  must  be  taken  to  have  adequate  spatial  resolution.  This 
fact,  combined  with  the  fact  that  collocation  methods  are  easy  to  apply 
in  the  nonlinear  case,  shows  that  pseudospectral  method  is  in  general 
preferable  to  Galerkin  or  Tau  methods. 
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TABLE  I. 


L2  Chebyshev  errors  for  the  solution  of  (4.6)  -  (4.9) 


N 

u  =  ( 1+x) u 

u  x 

Ut  =  (1'x)ux 

U  =  X  u 

t  X 

U  =  -X  u 
t  X 

17 

1.13  •  10_1 

vO 

1 

O 

• 

<r 

Ch 

1.16  •  io-1 

2.05  *  10-6 

33 

1.79  •  10'3 

4.7  •  10~7 

2.59  •  io“4 

1.05  *  10-7 

65 

8.5  •  10~5 

2.2  •  10"8 

1.22  •  io“5 

-9 

5  *10 
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